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We study the QCD phase diagram, in particular we study the critical points of the two main QCD 
phase transitions, confinement and chiral symmetry breaking. Confinement drives chiral symme- 
try breaking, and, due to the finite quark mass, at small density both transitions are a crossover, 
while they are a first or second order phase transition in large density. We study the QCD phase 
diagram with a quark potential model including both confinement and chiral symmetry. This for- 
malism, in the Coulomb gauge hamiltonian formalism of QCD, is presently the only one able to 
microscopically include both a quark-antiquark confining potential and a vacuum condensate of 
quark-antiquark pairs. Our order parameters are the Polyakov loop and the quark mass gap. The 
confining potential is extracted from the Lattice QCD data of the Bielefeld group. We scan the 
QCD phase diagram for different quark masses, in order to address how the quark masses affect 
the critical point location in the phase diagram. 
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Figure 1: A sketch of the QCD Phase Diagram, according to the collaboration CBM at FAIR It is 
usually assumed that the critical point for deconfinement coincides with the critical point for chiral symmetry 
restoration. 

1. Motivation 

Our main motivation is to contribute to understand the QCD phase diagram, for finite T and 
/I, to be studied at LHC, RHIC and FAIR. 

Using modem quark models for light quarks we now study chiral symmetry breaking, i. e. 
quark mass generation, at finite T. In the last Bormio Meeting we used the bottomonium and char- 
monium as good prototypes to study finite T quark-antiquark potentials. In that case it was suffi- 
cient to solve the Schrodinger equation with static lattice QCD potentials, since mh,mc » Aqcd 
and mh^nic » Tc allowed us to neglect in the quark sector, temperature effects, spontaneous chiral 
symmetry breaking, relativistic effects and coupled channels. However, to address light quarks at 
finite temperature T it is also necessary to include temperature, not only in the confining quark- 
antiquark interaction, but also in the quark mass generation and in the quark propagator. 

Thus the present work, not only addresses the QCD phase diagram, but it also constitutes the 
first step to allow us in the future to, 

- compute the spectrum of any hadron at finite T , 

- compute the interaction of any hadron-hadron system at finite T . 

Here we address the finite temperature string tension, the quark mass gap for a finite current 
quark mass and temperature, and the deconfinement and chiral restoration crossovers. We conclude 
on the separation of the critical point for for chiral symmetry restoration from the critical point for 
deconfinement. 



2. Fits for the finite T string tension from the Lattice QCD energy Fi 

At vanishing temperature T = 0, the confinement, modelled by a string, is dominant at mod- 
erate distances, 

V{r)^-^ + Vo + ar. (2.1) 
Izr 
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At short distances we have the Luscher or Nambu-Gotto Coulomb due to the string vibration + 
the OGE coulomb, however the Coulomb is not important for chiral symmetry breaking. At finite 
temperature the string tension cj(r) should also dominate chiral symmetry breaking, and thus one 
of our crucial steps here is the fit of the string tension o{T) from the Lattice QCD data of the 
Bielefeld Lattice QCD group, |g |3|, § |H]. 

The Polyakov loop is a gluonic path, closed in the imaginary time f4 (proportional to the 
inverse temperature T^^) direction in QCD discretized in a periodic boundary euclidian Lattice. It 
measures the free energy F of one or more static quarks, 

P{{))=Ne-^"l'^ , P''{0)P%r)=Ne-^'"^''^/'^ . (2.2) 

If we consider a single solitary quark in the universe, in the confining phase, his string will travel as 
far as needed to connect the quark to an antiquark, resulting in an infinite energy F. Thus the 1 quark 
Polyakov loop P is a frequently used order parameter for deconfinement. With the string tension 
a(r) extracted from the qq pair of Polyakov loops we can also estimate the 1 quark Polyakov loop 
P(0). At finite T, we use as thermodynamic potentials the free energy Fi and the internal energy 
U\, computed in Lattice QCD with the Polyakov loop |3|, ^ ^, ^. They are related to the static 
potential V{r) = —fdr with Fl{r) = —fdr — SdT adequate for isothermic transformations. In Fig. 
^ we extract the string tensions o{T) from the free energy F\{T) computed by the Bielefeld group, 
and we also include string tensions previously computed by the Bielefeld group [^]. 

We also find an ansatz for the string tension curve, among the order parameter curves of other 
physical systems related to confinement, i. e. in ferromagnetic materials, in the Ising model, in 
superconductors either in the BCS model or in the Ginzburg-Landau model, or in string models, 
to suggest ansatze for the string tension curve. We find that the order parameter curve that best 
fits our string tension curve is the spontaneous magnetization of a ferromagnet [^], solution of the 
algebraic equation, 

M (Tc M \ 

tanh — . (2.3) 



Msat \TM,a,/ 

In Fig. ^ we show the solution of Eq. 2.3 obtained with the fixed point expansion, and compare it 



with the string tensions computed from lattice QCD data. 
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Figure 3: Left: the mass gap m(0) solution of as a function of the quark current mass mo, in units of a = 1. 
Right: sketch of the effect of mo on the crossover versus phase transition of choral restoration at finite T . 



3. The mass gap equation with finite T and finite current quark mass mq. 



Now, the critical point occurs when the phase transition changes to a crossover, and the 
crossover in QCD is produced by the finite current quark mass mO, since it affects the order param- 
eters P or a, and the mass gap m(0) or the quark condensate {qq). The mass gap equation at the 
ladder/rainbow truncation of Coulomb Gauge QCD in equal time reads. 



m{p) 



o 

?M0 + — 7 
pi 



dk 



1 



271 ^B + m{kY 
pk pk 



pk 



pk 



{p-kf + ^l^ 



{p-kf + H^^ {p + kY + li^^ 2 ^{p + kf + ii^ 



{p + kf + H^ 
m{p)k 



m{k)p 



(3.1) 



The mass gap equation for the running mass m{p) is a non-linear integral equation with a 
nasty cancellation of Infrared divergences [I^, |ll]]. We devise a new method with a rational 
ansatz, and with relaxation, to get a maximum precision in the IR where the equation is nearly 
almost unstable. The solution m{p) is shown in Fig. |3|for a vanishing momentum p = 0. 

At finite T , one only has to change the string tension to the finite T string tension cj(r) [p^], 
and also to replace an integral in ft) by a sum in Matsubara Frequencies. Both are equivalent to 
a reduction in the string tension, a — a* and thus all we have to do is to solve the mass gap 
equation in units of o* . The results are depicted in Fig. ^. Thus at vanishing wiq we have a chiral 
symmetry phase transition, and at finite mo we have a crossover, that gets weaker and weaker when 
mo increases. This is also sketched in Fig. |[ 



4. Chiral symmetry and confinement crossovers with a finite current quark mass 

In what concerns confinement, the linear confining quark-antiquark potential saturates when 
the string breaks at the threshold for the creation of a quark-antiquark pair. Thus the free energy 
F (0) of a single static quark is not infinite, but is the energy of the string saturation, of the order of 
the mass of a meson i. e. of 2mo. For the Polyakov loop we get, 

P(0) ~A^£'"^"'»/^ . (4.1) 
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Figure 4: Sketches of the saturation of confinement (left), and of the corresponding crossover in the order 
parameter P polyakov loop (righ). 

Thus at infinite /mq we have a confining phase transition, while at finite mo we have a crossover, 
that gets weaker and weaker when mo decreases. This is sketched in Fig. ^. 

Since the finite current quark mass affects in opposite ways the crossover for confinement and 
the one for chiral symmetry, we conjecture that at finite T and ix there are not only one but two 
critical points (a point where a crossover separates from a phase transition). Since for the light u 
and d quarks the current mass mo is small, we expect the crossover for chiral symmetry restoration 
critical to be closer to the = vertical axis, and the crossover for deconfinement to go deeper 
into the finite /X region of the critical curve in the QCD phase diagram depicted in Fig. [I]. 
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